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1 Introduction 



A variety of methods, e.g. the exponential function [12], the hyperbolic tangent function Ii281 . 
the Jacobi elliptic function expansion 1261 . the sine-cosine function HI, the F-expansion ll24ll . the 
projective Riccati equation expansion HI, the ^-expansion Q, the ^)-expansion ll25l methods 
and their numerous extensions [2l|l|6l|9l[Ill[10l|2ll|27l[30l|3l][3l|Ml|3^ etc. to cite the main 
methods, are used to obtain analytical solutions for nonlinear differential equations. The common 
issue to all these methods is the rational form of the final solution. The rationality here is with 
respect to either trigonometric, hyperbolic, elliptic, exponential functions, or other given functions. 
Furthermore, the scope of these methods is restricted to polynomial differential equations. We 
describe in this paper the general principle from which all the above mentioned methods derive for 
the integration of autonomous nonlinear partial differential equations (NLPDEs). For the notations 
used in this work, see 1141 - flSl . 

2 Function expansion methods: general procedure 

Consider the constant coefficient partial differential equations 

f(^mW(x))=0, (2.1) 

where the nonzero positive integer s is the order of the equation and x= (x^ ,x") are independent 
variables. The dependent variable u = u{x) is a scalar valued function. If 

^("^'-"W) =-P("(*:,)[/^l]W,"(*:2)[^2]W,--- ,U{kr)[m^)) , (2-2) 

where r,ki,hi G N with max{/c,-, / = 1, • • • ,r} = 5, /j,- G {1,2, • • • and P is a polynomial whose 
the indeterminates are the r functions M(,(.,) [^^i] {^): the equation (12.11 ) becomes 

P {u{ko[f^i]{x),U(^k2)[h2]{x), - ■ ■ ,U{k,)[hr]{x)) =0 (2.3) 

which is called a polynomial autonomous partial differential equation. 

The change of variables u{x) =v (^) with ^ = aix^ +a2x'^-\ \-(Xnx", where a, G IR transforms 

(12.311 into an ordinary differential equation 

e(v(,,)[l](^),v(,,)[l](^),--- ,u^t,.)[m)) =0, (2.4) 

where Q is also a polynomial function whose the indeterminates are V(,(-,) [1] (^) ) ' = 1 ) " ■ " if. 
If the polynomial Q has only one monomial, then the equation (12.41 ) is simply solved by successive 
integrations. Thus, without loss of generality, one can assume that the polynomial 2 is a sum of 
at least two monomials. We propose in this paper to search for a solution of the equation (12.41 ) by 
using single, double or triple function expansion methods. 

2.1 Single function expansion method 

It consists to seek the function v = v(^), solution of the equation (I2.4I ). into the form 

v=A(F)+F(i)[l]B(F), (2.5) 
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where 

m m 

A{F) = £ aiF\ B{F) = £ b,F' 

i=—m i=—m 

and the function F = F (^) is a solution to an auxiliary differential equation which is in one of the 
forms 

F(i)[l]=7? or (F(i)[l])'=/?, (2.6) 

R being a rational function in F and m,/M G N, a,,Zj, G M. The higher order derivatives of F can be 
written into the form 

where the positive integer k>2 and R'^^, / = 0, 1, are also rational functions in F. 

One can adopt the following iterative process in practice for the determination of different pa- 
rameters of V. 



Estimation of the integers m, m. 

Let Ml ,M2, • • • ,Mv be the monomials of Q such that M\ contains the highest order derivative of the 
function v(^) and let M2 be a nonlinear or linear monomial. is linear only in the case where all 
remaining monomials are linear. Substitute in (12.41 ) the expression (12.51 ) of v(^) along with (12.61 ) and 
write the result in the form 

Q = \;Y.{Ki + F^,^[l]Li), (2.7) 
1=1 

where A',,L,, T are polynomials in F whose degrees depend linearly on m,m such that 

Mi = j{Ki + F^,)[l]Li). (2.8) 

Solve the system of linear algebraic equations obtained by balancing the degree of Ki with that of 
K2 and the degree of Li with that of L2 in F to determine the values of integers m,m. 

Estimation of the constants ai,bi, a,. 

Introduce into (12.71 ) the obtained values of m,m. Write the result in the form 

e = ^(^ + F(i)[l]L), (2.9) 

where K = ^4=1 Ki and L = ^4=1 ^i- Set to zero all coefficients of distinct monomials in K and L. 
This gives a system of algebraic equations whose the unknowns are the constants a,-, ft,- and a,. 

Remark 2.1. We have the tanh-expansion method if we take F{^) = tanh(^), the tan-expansion 
method if F(^) = tan(^), the exp-expansion method if F(^) = exp(^), the ^-expansion method if 

F{^) = where G"{^) = aG'{t,) + PG(^) with a, P G M. It is noticeable that in all these cases 
F'{t,) is a polynomial function in 7^(^) and hence a rational function as required. 



3 



2.2 Double function expansion method 

It aims at finding tlie function v = v(^), solution of the equation (12.41) . into tiie form 

v=A{F,G) + G(i)[l]C(f,G) +F(i)[l]G(i) [1]D(F,G), (2.10) 

where 

mi m2 m\ »i2 

A{F,G)= £ £ B{F,G)= £ £ 

i——m\ j——m2 i——m\ j——m2 

m\ m2 ni\ m2 

i——mi j——in2 /——mi j——in2 

and the functions F = F(^) and G = G(^) are solutions of an auxiliary system of differential equa- 
tions of the form 

F(i)[l]=/?i or (F(i)[l]) (2.11) 
G(i)[l]=/?2 or (G(i)[l])'=/?2, (2.12) 

Rk,k = \,2 being rational functions in F, G and mi,m2,mi,m2,mi,m2,'ni,m2 GN, Uij , j- , c,- , <i,- G 
M. The higher order derivatives of F and G can be written into the form 

F(,)[l] = /?f + F{i)[l]/?;'^ + G(i)[l]/?f + F(i)[l]G(i)[l]/?f, 
G(,)[l] = /?f + F(l)[l]/^^■*^ + G(l)[l]/^^''+F(l)[l]G(l)[l]4^ 

where the positive integer k>2 and 7?^*, / = 0, 1 , 2, 3, j = 1,2, are also rational functions in F, G. 

One can adopt the following iterative process in practice for the determination of different pa- 
rameters of V. 



Estimation of the integers m i , m2 , m i , m2 , m i , m2 1 ^ 1 1 ^2 • 

Let Mi,M2, ■ • • ,Mv be the monomials of Q such that M\ contains the highest order derivative of the 
function v(^) and let M2 be a nonlinear or linear monomial. M2 is linear only in the case where 
all remaining monomials are linear. Substitute in (12.41) the expression (12.101 ) of v(^) along with 
(I2.11I )- (I2.12I ) and write the result in the form 

2 = 7L(^'- + ^(i)[l]^^ + G(i)[l]5,- + F(i)[l]G(i)[l]7;), (2.13) 

where Ki , L, , 5, , , T are polynomials in F, G whose the degrees linearly depend on m 1 , m2 , m 1 , m2 , 
m[, m2, fhi, 1712 such that 

M,- = l(/:,- + F(i)[l]L,- + G(i)[l]5,- + F{i)[l]G(i)[l]7,-)- (2-14) 

Solve the system of linear algebraic equations obtained by balancing the degree of K\ with that of 
K2, the degree of Li with that of L2, the degree of with that of ^2 and the degree of Ji with that 
of J2 in F and G to determine the values of integers m\,m2,mi,m2,ini,rn2,mi,rh2- 
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Estimation of the constants aij,bij,Cij,dij, a, . 

Introduce into (12.131) the obtained values of m\,m2,rhi,m2,rh\,m2,ni\,fh2- Write the result in the 
form 

e = ^(^ + F(i)[l]L + G(i)[l]5 + F(i)[l]G(i)[l]7), (2.15) 

where K = LjLi^i) L = Y^i^^Li, S = L/Li'^r and J = LLi-^i- Set to zero all coefficients of 
distinct monomials in K, L, S and /. This gives a system of algebraic equations whose the unknowns 
are the constants aij,bij,Cij,dij and a,-. 

Remark 2.2. We have the (sinh, cosh) -expansion method if we take F(^) = sinh(^) and G(^) = 
cosh(^), the (sin, cos) -expansion method if F(^) = sin(^) and G(^) = cos(^). We obtain a more 

general formulation of the -expansion method by setting F(^) = and G(^) = jjj^, 

where H"(^) + XH{i,) = /u with X,^ G M, and a more general formulation of the projective Riccati 
equation expansion method if we choose F(^),G(^) such that F'(^) = aF(^)G(^) and G'(^) = 
lj + a^G^(t,) - pF(^) with a,p,^ G M^ In all these cases, F'{t,) and G'(^) are also polynomial 
functions in F(^),G(^) and hence rational functions as required. 

2.3 Triple function expansion method 

Here, we need to express the function v = v(^), solution of the equation (12.41 ). into the form 

V = A + b1'1F(i)[1]+b2=1G(i)[1]+b3'^i)[1]+c1-'F(i)[1]G(i)[1] 

+ c2'1F(i)[1]//(i)[1]+c3''G(i)[1]//(i)[1]+DF(i)[1]G(i)[1]//(i)[1] (2.16) 

where A,B^\B^\B^\C^\C^\C^\D are {unctions of F,G,H givenhy 

mi m2 1113 '"11 ""l-Z '"1.3 

^= L L L ^u^p'G^fi\ Bi'i= III b\jjj^'Gm\ 

[——mi j——m2k——m^ — '"Li j——'^l2^——'^i.3 

1712 1 ^2 2 W^2 3 '^3 I '^^3 2 '"3 3 

B''= LIE 5^'^= I x x 

i——m2A j——^l,2 k——rn2^3 '——"^3.1 7— — "^3,2 ^— — W3.3 

mi 1 mi 2 3 "'2 1 2 3 

c^'^= III c2'^= III cy.i^w, 

i=—m\ I j=—rh\ 2k=—m\ 3 '«2.l j=^'fi2.2k=— 1112,3 

'«3.l '"3.2 'f'3,3 jfii mo m3 

C''= III ^:hF'GJH\ 0= I I I 

i=— ms.l j=— 11132 k=—fh3,3 i=—m\ j=—m2 k=—m3 

and mi,mi,mij,mij are positive integers, aij\ic,b\'^j i^,cf^j i^,dij,k G K; the functions F = F(^), G = 
G(^) and H = H{i^) are solutions of the following auxiliary system of differential equations: 

F(i)[l]=/?i or (f(i)[l])2=/?i, (2.17) 
G(i)[l]=/?2 or (G(i)[l])'=/?2, (2.18) 
//(i)[l]=/?3 or (//(i)[l])'=/?3, (2.19) 
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Rk, k = 1,2,3 being rational functions in F, G, H. The higher order derivatives of F, G and H can 
be written into the form 

f(,)[l] = <^ + F(i)[l]/?;'^ + G(i)[l]/?f + //(i)[l]/?^'^ + F{i)[l]G(i)[l]/?f 

+ F(i)[l]//(i)[l]/?f + G(i)[l]//(i)[l]/?f + f(l)[l]G(l)[l]//(l)[l]/^]•^ 

G(,)[l] = <^ + F(i)[l]4^ + G(i)[l]/?f + //(i)[l]/?^'^ + F(i)[l]G(i)[l]4^ 

+ F(i)[l]//(i)[l]/?^'*^ + G(i)[l]//(i)[l]/?f + f(l)[l]G(l)[l]//(l)[l]/^^•^ 

= <' + F(i)[l]4*^ + G(i)[l]/?f + //(i)[l]/?^'*^ + F(i)[l]G(i)[l]4*^ 

+ F(i)[l]//(i)[l]/?f + G(i)[l]//(i)[l]/?f + F(i)[l]G(i)[l]//(i)[l]/?^'\ 

where the positive integer k>2 and R''^ , i = 0, 1,2,3,4,5,6,7, j = 1,2,3, are also rational func- 
tions in F, G, H. In practice, the following iterative process leads to the determination of different 
parameters of v. 

Estimation of the integer mi.mi.rhi j.nii j. 

Let Ml ,M2, • • • ,Mv be the monomials of Q such that M\ contains the highest order derivative of the 
function v(^) and let M2 be a nonlinear or linear monomial. M2 is linear only in the case where 
all remaining monomials are lineai\ Substitute in (12.41) the expression (12.161 ) of v(^) along with 
(12. 17I )- (I2.19I ) and write the result in the form 

Q = 7 L + ^(1) [1] l^U + [l]L2,i + //(I) [1] L3,,- + [l]G(i) [1] 
1=1 

+ T^i(i^(i)[l]//(i)[l]% + G(i)[l]//(i)[l]53,,- + F(i)[l]G(i)[l]//(i)[l]7,-), (2.20) 
1=1 

where Kt, Ljj, Sjj, Ji, T are polynomials in F,G,H whose the degrees linearly depend on m,, m,-, 
rhij, ihij and 

Mi = ^ (/:,- + F(i)[l]Li,,- + G(i)[l]L2,,- + //(i)[l]L3,,- + F(i)[l]G(i)[l]5i,-) 

+ i(F(i)[l]//(i)[l]52,,- + G(i)[l]//(i)[l]53,,- + F(i)[l]G(i)[l]//(i)[l]70^ (2.21) 

Solve the system of linear algebraic equations obtained by balancing the degree of Ki with that of 
K2, the degree of Ljj with that of Ly 2, the degree of Sj^i with that of Sj^i and the degree of Ji with 
that of J2 in F, G and H to determine the values of the integers mi,mi,mi j,mi j. 

Estimation of the constants aij^k, c-^j dij^k, o^r • 

Introduce into (12.201 ) the obtained values of mi,mi,mij,mij. Write the result in the form 

Q = l(/: + F(i)[l]Li+G(i)[l]L2 + //(i)[l]L3+F(i)[l]G(i)[l]5i) 

+ ^ [l]//(i) [1] 52 + G(i) [l]//(i) [1] 53 + [l]G(i) [l]//(i) [1] J) . (2.22) 
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where ^ = I)Li^;, Lj = i:i^,Ljj, 5y = Ij^i 5,-; and 7 = H^i 

Set to zero all coefficients of distinct monomials in K, Lj, Sj and /. This gives a system of algebraic 
equations whose the unknowns are the constants at^j^kMi] k^^Yj k^^Uj^k ^'^'^ ^i- 

Remark 2.3. We have the (i'?i,c«,(i?i) -expansion method if we take F{^) = sn(^), G(^) = cn{^) 
and H(^) = dn(^), where sn(^) = sn(h„k), c«(^) = cn(h„k) and dn(i^) = dn{^,k) with < ^ < 1 
are the basis Jacobi elliptic functions [7 |. The function sn, cn and dn are solutions of the first order 
ordinary differential equations 

(w'(^))' = (l-w2(^))(l-^V(^)), (2.23) 



w 



'(^)) =(l-w2(^))(^V(^) + l-^2), (2.24) 



(w'(^))' = (1 + {w\^)+k^ - 1) , (2.25) 

respectively. When ^ — )■ the Jacobi functions degenerate to the trigonometric functions, i.e., 

Iim5?i(^,^) = sin(^), limc?i(^,^) = cos(^), lim<i?i(^,^) = 1. 

k^O k^O k-^Q 

When — )• 1, the Jacobi functions degenerate to the hyperbolic functions, i.e., 
Urn sn{l„ k) =tanh{t,), limcn(^,fc) = — ttfT) limdn{t„k) 



k^i yt-s-i cosh(^)' k^i ^ cosh(^)' 

The function sn, cn and dn have the algebraic properties 

sn'^{t„k)+cn^{^,k) = 1, k^sn^{^,k) + dn^{^,k) = 1, 

k^cn^{l„k) + \-k^ = dn^{i„k), cn^{l,,k) + {l-k^) sn^(X,k) = dn^{^,k) 
and the differential properties 

--^sn{^,k) = cn{^,k)dn{^,k), --^cn{^,k) = —sn{^,k)dn{^,k), 
■^dn{^,k) = —k^ sn{^,k) cn(^,k) . 

3 Application to some relevant NLPDEs in physics 

(A) Case 1: Burger-Fisher equation 

Let us use the single function expansion method to analyze the Burger-Fisher equation 132 



Uxx + UUx — Ut + U — U^ = 0, (3.1) 

where u = u{t,x). We first substitute into (13.11 ) the variables u{t,x) = v(^), ^ = + (3x, with 
a, (3 G M to obtain 

P^v^^ + Pvv^-av^ + v-v^ = 0. (3.2) 
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We take the expansion (12.51) to look for the solution of the equation (13.21) with the assumption 
that the function F = F{b,) satisfies the auxiliary equation f '(^) = 1 — f ^(^). Balancing the 
degree of the nonlinear term vv^ with that of the linear term v^^ yields m = m = l. This allows 
to take the ansatze 

V© = £ atF'{i,)+F'g) t bjF\i,), (3.3) 

i=-\ j=-\ 

where ai^bj are constants to determine. Here F(^) = tanh(^) and (13.31 ) becomes 



v(^) = -bi tanh^ (^) - bo tanh2(^) + {ai+bi-b-i) tanh(^) + ao + - 



tanh(^) 



(3.4) 



Substitute (13.41) into (13.21) . Then, collect and set to zero the coefficients of each power 
of tanh'(^). Solving the resulting system with respect to the parameters a;,Zjy,a,P, several 
solution sets are obtained leading to the following solutions for the equation (13.11 ): 8 G 



Ml [t,X 

UT,{t,X 
U4{t,X 
Us{t,X 

U(i{t,X 
UT{t,X 
U^{t,X 



0, 

1, 
1 

- + 



2 tanh(|f)' 

1 £ /£ 
- + - tanh -t 

2 2 \2 



1 



+ ■ 



2 tanh[|(ff+x)] 



1 £ u 

- + - tanh 

2 2 



4 V2 



-t+x 



i + |tanh(|n + 



tanh (10' 



1 £ u 

- + - tanh 

2 4 



£ /5 



8 V 2 



-t+x 



+ 



tanh [| (ff+x)] ' 



(B) Case 2: Burger-Fisher equation 

Consider now the double function expansion method for the analysis of the Burger-Fisher 
equation B21 

Uxx + uUx — Ut + u — u'^ = 0, (3.5) 

where u = u{t,x). We first substitute into (13.51 ) the variables u{t,x) = v(^), ^ = + Px, with 
a, P G M to obtain 

(3^v^^ + (3vv^-av^ + v-v^ = 0. (3.6) 

We take the expansion (12.101 ) to look for the solution of the equation (13.61 ) assuming that the 
functions F = F(^) and G = G(^) satisfy the auxihary equation F'(^) = G(^) and G'(^) = 
F(^). Balancing the degree of the nonlinear term vv^ with that of the linear term v^^ yields 



8 



mi = m2 = rhi = m2 = rhi = m2 = mi = m2 = 2. This allows to take the ansatze 



^=11 {a/j + ^-j^(i)[l]+c,-jG(i)[l]+J,yF(i)[l]G(i)[l]}f'G^ 

i=-2j=-2 



(3.7) 



where ai j,bij,Cij,dij are constants to determine. Note that F(^) = cosh(^) and G(^) = 
sinh(^). Substitute (13.71) into (13. 6I ). collect and set to zero the coefficients of each power of 
cosh'(^) sinh^(^). Solving the resulting system with respect to the parameters aij, bij, Cjj, 
dij, a, (3, several solution sets are obtained leading, in addition to the solutions found in the 
Case 1, to the following solutions of equation (13.51) . (a is an arbitrary constant): 

ug{t,x) 
uio{t,x) 



a [sinh(2f + cosh(2f , 
1 +a [sinh(f +x) +cosh(f +x)] , 



un{t,x) 
uu{t,x) 



1 



sinh [ t + -x\ + cosh \t + -x 



1 



l+a 



n 2 



sinh ( + j + cosh f + ^-^ 



(C) Case 3: Fifth-order KdV equations 

Consider the well known fifth-order KdV (fKdV) equations ||32]| in its standard form: 

Ut + au^Ux + 5UxU2x + PUU3X + "Sjc = 0, (3.8) 

where a,5,p are arbitrary nonzero real parameters and u = u{t,x) is a sufficiently smooth 
function. A variety of the fKdV equations can be retrieved from this equation by changing 
the real values of the parameters a, 5 and p. However, five well known forms of the fKdV 
equations are of particular interest in the Uterature. There are: 

(CI) The Sawada-Kotera (SK) equation ||29l given by 

Ut + Su^Ux + 5uxU2x + iuuj,x + M5jc = 0; (3.9) 
(C2) The Caudrey-Dodd-Gibbon equation ||8l given by 

Ut + ISOm^m.v + 30uxU2x + 'iOuu^x + u^x = 0; (3.10) 
(C3) The Lax equation |[23l provided by 

Ut + 30u^Ux + 20uxU2x + lOuu^x + U5x = 0; (3.11) 
(C4) The Kaup-Kupersmidt (KK) equation ||22l expressed as 

Ut + 20u Ux + 25 UxU2x + ^Ouu^x + usx = 0. (3.12) 
(C5) The Ito equation |fT9l written as 

Ut +2u^Ux + 6UxU2x + 3llU3x + ll5x = 0. (3.13) 
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avv + p^vsy + p^vvsy + ^ov\ + p^6vj,v2y = (3. 16) 



It is important to note that there is another significant fifth-order equation that appears in the 
literature in the form 

Uf +OUUx + 5ut,x — pu5x = 0, (3.14) 

where a, 5, p are constants. This equation is called the Kawahara equation. The standard form 
of the Kawahara equation [|20,I is a fifth order KdV equation of the form 

Ut + 6uUx + ut,x — = 0, (3.15) 

that describes a model for plasma waves, capillary-gravity water waves. The Kawahara equa- 
tion appears in the theory of shallow water with surface tension and in the theory of magneto 
acoustic waves in a cold collision free plasma. In [ 13J we have studied the Kawahara equation 
(13.141 ) using the symmetry group analysis method. 

In the new variables y = at + and v{y) = u{t,x), (13.81 ) is reduced to the ordinary differen- 
tial equation 

fv,, _lR5,,. _lR3„,„,_ _lR^,,2,, _lR3x-^^^ 

while equation (13.141 ) is reduced to 

avy + pavvj, + p^Svay - p^pv5>. = 0. (3. 17) 

The single function expansion method suggests to take the ansatze 

v{y) = £ ar{y)+F'{y) £ bjFJ{y) (3.18) 

i=-2 j=~2 

for the equation (13.161 ) and the ansatze 

4 4 

v{y) = £ aiF\y)+F'iy) £ bjF^iy) (3.19) 

r=-4 j=-4 

for the equation (13.171 ). where ai,bj are constants to determine and F is a sufficiently smooth 
function. Here, we assume that the function F = F{y) satisfies the auxiliary equation F'{y) = 
l-F^{y), i.e. F(y) = tanh(3;). Thus, (l3?T8l) becomes 

v(3;) = -b2 tanh'^(3;) - bi tanh^ (y) + (a2 - + ^2) tanh^ (y) + {ai - b^i +bi) tanh(y) 

+ '-T^ + '-^ (3-20) 
tanh(3;) tanh^(3;) 

and (13.191 ) becomes 

v{y) = -Zj4tanh^(3;) -^3tanh^(3;) + (^4 -Z?2 + ^4)tanh'*(3;) + {as-bi +b3)tmh^{y) 
+ (a2-^() + ^2)tanh^(y)+ao-^-2 + ^o + -""^ ^-3+^-1 



tanh(j) 

tanh^(y) tanh''(y) tanh'*(j) 

Substitute (13.201 ) into (13.161 ) and (13.211) into (13.171 ). Then, collect and set to zero the coeffi- 
cients of each power of tanh(y). Solving the resulting system of algebraic equations, several 
sets of parameters {ai,bj,a,^} are obtained leading to the following results: 
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(i) Analytical solutions to the Kawahara equation (13.151 ): 



Ml {t,x) 
U2{t,x) 

Ui,{t,x) 

U4{t,x) 



ao', 



ao 



ao 



ao + 



35 



35 
169 
35 
169 
35 



tanh^(cp(f,x)) - -tanh^(cp(f,;c)) 



1 



1 



1 



5408 



tanh^(cp(?,x)) 2 tanh**(cp(?,x)) 
tanh'^(\|/(f,x)) + ^ 



where (p{t,x) 



1352 

3\/T3(338qo-23)^ 



tanh^(\|/(f,x)) + 



tanh'*(\|/(?,x)) 
1 



tanh^{\\f{t,x)) 



/13 



X and 



4394 ' 26 

(ii) Analytical solutions to the Sawada-Kotera equation (13.91) : 



3^/n(2704ao-9) ^ 
70304 ' 52 



U2{t,x) 
Ui,{t,x) 

M4(f,x) 
M5(?,x) 

U(,{t,x) 

uj(t,x) 



ao; 



8p2-i2(32tanh2(q)(f,x)); 



8p2- 12p2 
ao - 6p2 



1 



tanh^((p(?,x)) ' 

ao-6P^tanh^ (cp(f,x)); 
1 



tanh^ ((p(f,x)) ' 



8p2-12p2 



tanh^(\|/(?,x)) + 



1 



ao - 6p2 



tanh^ + ' 



tanh^(\|/(f,x)) 
1 



tanh^ (\|/(f,x)) 



where 



cp(f,x) = 16p^f - Px, \|/(?,x) = 256p^? - px, 
(p(f ,x) = p (76p^ + ag - 40p2ao) t - px, \^(?,x) = p ( 16p^ + 5a,^ - 40p2ao) 
(iii) Analytical solutions to the Caudrey-Dodd-Gibbon equation (13.101 ): 



Ul{t,X 
U2{t,X 

U3{t,X 

U4{t,X 
Us{t,X 

U(i{t,X 

Uj{t,X 



ao; 



-p2-2p2tanh2((p(?,x)); 

lp^-2p^^i ; 

3 tanh^(cp(?,x)) 

ao-P^tanh^ (\|/(f,x)) ; 
«o-P' 



1 



tanh^ (\|/(f,x)) ' 



-p2-2p2 



tanh^ {(^{t,x)) + 



1 



«o-p' 



tanh"^ {\\r{t,x)) + 



tanh^ {(p{t,x)) 
1 



tanh (\|/(f,x)) 
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where 



(?{t,x) = 16p^f - ^x, (^{t,x) = 256p^f - ^x, 
y\f{t,x) =4(3 (l9pV45ag - 60fao)t-^x, =4p (4pV45af) - eOp^ao) t - ^x. 

(iv) Analytical solutions to the Lax equation (13.11b : 



Ml {t,X 
U2{t,X 

UT,{t,X 

U4{t,X 
U5{t,X 

U(,{t,X 

My {t,X 

U^{t,X 

U(){t,X 



ao', 



4p2-6p2tanh2((p(f,x)); 



4p2_6p2 



1 



tanh^(cp(f,x)) ' 
ao - 2P^tanh^ ; 
1 



4p2_2p2 
4p2 - 2p2 
4p2_6p2 



tanh^ ' 
tanh^ {^{t,x)) + 3 



1 



3tanh^ ((])(?, x)) + 
tanh^ {(^{t,x)) + 
tanh^ (^(^■^)) + 



tanh^ {(^{t,x))_ 
1 

tanh^ ((])(?, x)) 
1 

tanh^ (cp(f,x)) 
1 



tanh^ 



where 



(p(f,x) = 56p^f-px, (p(?,x) = 896p^f-px, ^(f,x) = 336p^f - px, 
,x) = 2p (28p4 + 15ag - AOfao) t - px, \i7(?,x) = 2p (48p4 + 15flg - AO^'^qq) t - px. 
(v) Analytical solutions to the Kaup-Kupersmidt equation (13.121 ): 



Ul{t,X 
U2{t,X 

UT,{t,X 

U4{t,X 
U5{t,X 

U(,{t^X 

Ui{t,X 



^^-^^heinhHi?it,x)y, 

p2-lp2 i ; 

2 tanh^(cp(f,x)) ' 
8p2- np^tanh^ (\|/(f,x)); 

1 



8p2 - 12p2 



tanh^ (\|/(f,x)) ' 
tanh^ (cp(f,x)) + 



1 



tanh (cp(f,x)) 



8p2 - 12p2 



tanh^ {^{t,x)) + ■ 



1 



tanh (\|/(f,x)) 



where 



(p(?,x) = p^f-px, cp(f,x) = 16p^?-Px, \|/(;,x) = 176P^r-px, \|/(f,x) = 2816P^?-px. 



12 



(vi) Analytical solutions to the Ito equation (13.13b : 



Ml 


{t,x) 


= ao; 




[l,X) 


— tp 


M3 


it x] 


= 4B^- 


U4 


{t,x) 


= 20p2- 


Us 


{t,x) 


= 20p2- 


"6 


t,x) 


= 4p2- 


M7 


{t,x) 


= 20p2- 



1 



tanh2(px)' 
)p2tanh2 (96p5f 
1 



tanh^ (96p5f - px) ' 
1 



tanh"^ (P^) + ■ 



tanh^ (Px; 
tanh2(l536p^f-px) + 



tanh'^(1536p5;-px) 



Finally, let us emphasize that all these methods involve cumbersome computations and hence 
do require the use of powerful computers. Except for this disadvantage, our analysis can be straight- 
forwardly extended to multiple function expansion methods and to the investigation of systems of 
partial differential equations. 
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